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It is well-known that allowing for spatial curvature affects constraints on cosmological parameters 
such as the dark energy equation of state parameters. Here we study the effect of curvature on 
constraints on parameters used to test general relativity (GR) at cosmological scales, commonly 
known as modified growth (MG) parameters, as while current data taken in the context of the 
ACDM model points to a universe that is flat or very close to it, this constraint may not hold in 
modified theories of gravity. Using the latest cosmological data sets we find that MG parameters 
are correlated with the curvature parameter Qk and the constraints on the MG parameters are 
weakened compared to when n& is not included in the parameter analysis. We next use various 
, future simulated data sets, including cosmic microwave background, weak lensing, and Integrated 

Sachs- Wolfe-galaxy cross-correlations, where the fiducial model is spatially curved but we assume 
a flat model when fitting the MG parameters. We find the assumption of a spatially flat model 

Qon a spatially curved universe does indeed cause an artificial shift in the constraints on the MG 
parameters, in some cases even producing an apparent deviation from GR in the MG parameter 
f"*; ' space. For our simulated data, tension with GR begins to manifest itself for fiducial models with 

| > 0.02 and apparent deviations appear for \Qk\ > 0.05. We find that for negatively curved 

models the apparent deviation is more significant. The manifestation of this apparent deviation from 
i i GR due to the assumption of spatial flatness above leads one to conclude that, when using future 

high-precision data to perform these tests, spatial curvature must be included in the parameter 
, analysis along with the other core cosmological parameters and the MG parameters. 

^ 1 PACS numbers: 95.36.+x,98.80.Es,98.62.Sb 

P 1 I. INTRODUCTION 

ryj ■ 

In an era of precision cosmology several studies have shown that allowing curvature can affect constraints on other 
cosmological parameters, particularly the equation of state of dark energy (see for example Refs. [U-IH). With the 
{SJ ' current interest in constraining modified growth (MG) parameters to test general relativity on cosmological scales, 
. one might wonder what effect curvature will have on these constraints and their implications. 

Indeed the effort of understanding the cause of cosmic acceleration by measuring the equation of state of dark 
energy has been joined by studies in constraining deviations from general relativity at cosmological scales [7H32| and 
most recently by using current data to constrain the MG parameters [33U4lj to determine if cosmic acceleration could 
be due to some extension to general relativity rather than an unknown dark energy component pervading the universe. 
So far, all constraints on MG parameters have shown that current observations are consistent with general relativity. 
However, these current constraints are not stringent enough to provide any definitive answer. With the exception of 
Ref. [36j . to date these studies have been carried out assuming a fiat universe. It is worth exploring what impact, if 
any, allowing curvature has on both constraining the MG parameters and the performance of these tests. 

In this paper, we explore the effect of spatial curvature on testing general relativity (GR) at cosmological scales. We 
do this because though current observations when interpreted via a Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe governed by general relativity point to a universe that is flat or very close to it, modified theories of gravity 
& • or inhomogeneous cosmological models may in fact require curvature to fit observations [42H44} . In order to perform 
the analysis, we expand the framework of modified growth (MG) equations and parameters. We also expand the 
corresponding numerical framework described in [Ibl ] to include curved models and apply these changes to the publicly 
available code, Integrated Software in Testing General Relativity (iSiTGR) 45]. To constrain the MG parameters 
while allowing for spatial curvature, we will then use current observations including: the WMAP7 temperature and 
polarization spectra (CMB) (46[; the matter power spectrum (MPK) from the Sloan Digital Sky Survey (SDSS) DR7 
[47j : the Integrated Sachs- Wolfe (ISW)-galaxy cross-correlations 0, SH; the refined HST-COSMOS weak- lensing 
tomography [50j |; WiggleZ BAO measurements [5l|; the supernovae compilation of the Supernovae Cosmology Project 
(SCP) 1531 and references of other compiled supernovae (SN) therein; the prior on Hq = 74.2 ± 3.6 km/s/Mpc given 
by Ref. (53|; and finally, a prior on the age of the Universe (AGE) 10 Gyrs<AGE< 20 Gyrs. We will explore the effect 
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this has on the constraints compared to the flat case. Next, we will analyze the correlations between the curvature 
parameter f2^ and the MG parameters. Finally, we will explore what effect assuming a flat model when curvature is 
present will have on the best-fit MG parameters, to see if a bias is introduced by such an assumption and discuss how 
this bias is related to the correlation coefficients. We conclude in the last section. 



II. GROWTH EQUATIONS INCLUDING SPATIAL CURVATURE IN GENERAL RELATIVITY 

The perturbed FLKW metric written in the general conformal Newtonian gauge is given by 

ds 2 = a(r) 2 [-(l + 2^)dr 2 + (1 - 2(p)^ ij dx i dx j ], (1) 

where </> and "0 are scalar potentials describing the scalar mode of the metric perturbations, r is conformal time, a(r) 
is the scale factor normalized to one today, and the Xj's are the comoving coordinates. 7y- is the 3- metric, which can 
be written in coordinates (x, y, z) as 



7y = 



k r 2 

l + — (x 2 + y 2 + z 2 ) 



(2) 



where K — —Q,i-T-L 2 is the spatial curvature, and ~Hq is the Hubble parameter today. 

As discussed in pH . l55j . when working in a non-flat universe the Fourier modes are generalized as eigenfunctions, 
G, of the Laplacian operator such that 

V 2 G(k,x) = -k 2 G(k,x). (3) 

In our analysis we expand perturbations in terms of G and its spatial covariant derivatives (denoted by | ) as seen in 
d 

Two very useful equations can be obtained using the first-order perturbed Einstein equations. These equations 
relate the scalar potentials to one another, as well as both the gauge- invariant, rest-frame overdensity, A,, and the 
shear stress, Oi (where i denotes a particular matter species). First combining the time-space and time-time equations 
gives the Poisson equation, while the traceless, space-space component of the equations gives a relation between the 
two metric potentials. Explicitly, these equations are 

(fc 2 - 3K) <(, = -AirGa 2 ^ Pl A t (4) 

i 

fc 2 (^-0) = -l2^Ga 2 Y J P^^ + w^)^ (5) 

i 

where pi is the density for matter species i. 

The quantity A^ is important because its evolution can be used to describe the growth of inhomogeneities. It is 
defined in terms of the fractional overdensity, Si = Spi/p, and the heat flux, qi (g, is related to the divergence of the 
peculiar velocity, 9i, by 9i = ), as 

Ai=8i + 3n^, (6) 
k 

where % = a/a is the Hubble factor in conformal time. By considering mass-ener gy c onservation, the evolution of A 
for an uncoupled fluid species or its mass average for all the fluids is given by (4p| |56|| 

,2 „-#M 1 



A = 3(1 + w) \J> + Hil)) + 3HwA - \ k z + 3 [H z - nj j ~ - 3H(1 + w)a, (7) 

where w — p/p is the equation of state of the fluid. Combining Eqs. (HJ, <j5j) , and along with the evolution 
equations for a(r), the growth history of structures in the universe can be fully described. 



III. MODIFICATIONS TO THE GROWTH EQUATIONS FOR SPATIALLY CURVED MODELS 
A. Modified growth equations in the conformal Newtonian gauge 



Recently, one of the major routes to testing general relativity has been parametrizing both the relation between 
the two metric potentials <j) an d ip m the perturbed FLRW metric (an inequality in this relation has been called 
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gravitational slip by J17l) as well as modifications to Poisson's equation, Eq. (|4|. Examples of this approach can 



be seen in to name a few. To date though, explorations of these modifications have only been 

performed in flat spacetimes. We now focus on extending the modified growth equations seen in, for example, [33|, H3] 
to non-flat cases. 

Extending the parametrized modifications of the growth equations, ((4]) and ([5]), proposed by 33] to non-flat models 
gives 

(k 2 -3K)4> = -AirGa 2 J2p i A l Q (8) 

i 

k 2 {iP-R(j)) = -U-KGat^piil+WijaiQ, (9) 

i 

where Q and R are the modified growth parameters (MG parameters). A modification to the Poisson equation is 
quantified by the parameter Q, while the gravitational slip is quantified by the parameter R (at late times, when 
anisotropic stress is negligible, tp = R<P)- As discussed in our earlier paper [Icj we use the parameter T> = Q(l + R)/2 
instead of R not only to avoid a strong degeneracy between Q and R, but also to have a parameter which can be 
directly probed by observations. To obtain a modified growth equation written in terms of only Q and T>, we combine 
Eqs. ||SJ} and ©, giving 

fc 2 (V + </>) = ^g^ E^ -127rGfl 2 ^> t (l + Wl KQ- (10) 



B. Modified equations in the synchronous gauge 

In order to perform the tests, we must implement the MG framework into numerical codes that allow comparisons to 
the data and calculations of parameter constraints. This is done by an extended version, with the inclusion of spatial 
curvature, of the publicly available package ISiTGR which is an integrated set of modified modules for the publicly 
available codes CosmoMC [(32| and CAMB [63] ■ CAMB is used to calculate the various CMB anisotropy spectra [Cj T , 
Cj E , Cf , Cf B ) as well as the three-dimensional matter power spectrum Ps(k, z) all of which are very powerful in 
constraining both the growth history of structure in the universe as well as the expansion history of the universe, and 
thus MG parameters. 

The package CAMB is written in the synchronous gauge, where the perturbed FLRW metric is written as: 

ds 2 = a(T) 2 [-dT 2 + ( 7ij + h^dx'dx ]. (11) 

Thus, instead of using the metric potentials <j> and tp of the conformal Newtonian gauge, it uses the metric potentials h 
and r\ consistent with the notation of 55] . These new perturbation quantities are related to the trace (h) and traceless 
(h + 6r/) part of the total metric perturbation hij. Explicitly, by writing hij expanded in terms of G, described above, 
for a single mode we have [HI] 

hj = ~HjG +(h + 6 V )(k- 2 G {lJ + pijG). (12) 

Now combining the perturbed Einstein's equations in the same way as discussed to get Eqs. (j4]) and <j5j> , we have for 
the synchronous gauge [64| 

(k 2 - 3K) (r; - Ua) = -47rGa 2 ^ / 9 l A J , (13) 

i 

k 2 (a + 2Ha--q) = -12nGa 2 ^ p t {l + w l )cr l , (14) 

i 

where a = (h + 677) /2fc 2 . 

Armed with the knowledge that both Aj and Ui are invariant between these two gauges and quickly comparing Eqs. 
(|3]) and (fHij) as well as ([5]) and {HP , we can see the metric potentials in the two gauges are related to one another by: 



4> = i] — "Ha, 
tp = a + Ha, 



(15) 
(16) 
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CAMB evolves the metric potential rj (this is done by evolving the quantity krj). For each matter species, the code 
also evolves: the matter perturbations, Sf, the heat flux, and the shear stress <7j according to the synchronous 
gauge evolution equations given in }56j |. The code has been written in such a way that the evolution of all other 
variables is changed simply by adjusting the evolution of the metric potential r\. So the first step to modifying CAMB 
to be consistent with the modified growth Eqs. and (H~0|) is to derive an equation for the evolution of 77 consistent 
with those equations. We begin by subbing (p~5)) into ((5} and taking the time derivative. This gives 



k 2 f ] = — —Y 

1 2K n ^ 



Pi(a)(QAi + QAi) + QAi—pi(a) 

dr 



Hk 2 a + Hk 2 a, 



(17) 



where we have used the condensed notation, Pi(a) — RirGa 2 pi and Kf 1 = 1 — 3K/k 2 . Next an expression for a is 
obtained by substituting Eqs. © and (jTHJ) into Eq. (fit))) , yielding 



a = -Ha-^^pi{a 



2V-Q 



/1 



3Q(1 + Wi)<7i 



(18) 



Now subbing the time derivative of pi (obtained from matter conservation) as well as Eqs. ([7]) and (fT8|) into (Ti~7]) 
gives: 



k 2 rj = 



-J-^^aJ QA l -HQA t + 3Q(l + w l ) (i + Hjf) -3^(1 + ^)^ 
-Q/i| + 2HVA, - UQAi + 3HQ K fl (l + tufa \ - {U 2 - U)k 2 a, 



where 

fi = k 2 +3(H 2 -%). 

Next we substitute Eq. (fl~6|) and the time derivative of Eq. (TTSl) into Eq. (fT9| , giving 
k 2 i) = 



(19) 



(20) 



Q^i - UQAi + 3Q(1 + w t ) (r) + (H 2 - H)aj - 3H(1 + w l )a t 
-Qhj + 2HVA i -HQA i + mQK fl {l + Wi)(Ti \ - {U 2 -U)k 2 a. 



(21) 



Since the variables in CAMB are evolved in the synchronous gauge, we must make sure that all the quantities in the 
■q are the synchronous gauge quantities. We know that A and a are gauge invariant, so we need not be concerned 
about them. However, qi in the above equations is still a conformal Newtonian gauge-quantity. Using Eq. 27b of [56[ 
(first converting from 8 to q as described above) the transformation of qi between these two gauges is written as 



q^ = 4 B) + (1 + m)ka. 

Upon subbing Eq. (|22[) into Eq. (1211) and collecting the terms we finally arrive at an equation for -q: 

2U[V-Q} + Q) Ai - Q(l + Wl )k 2 a 



(22) 



V 



2K 



-Qfif- + 3HQ(K fl - 1)(1 + wt)*, \ - 



n 2 -u 

fQ 



k 2 a, 



with: 



f Q = k 2 



3Q 
2K fl 



(23) 



(24) 



5 



Now as discussed in the time derivative of the sum of the Newtonian metric potentials, <j>+ip, goes into evaluating 
the contribution of the ISW effect in the CMB temperature anisotropy. This quantity needs to be consistent with the 
modified growth equation (flQ|) . Thus, simply taking the time derivative of ([TO]) and subbing in for A and pi gives: 



97? 

{l+3wi)Q + — )H-Q 



3(1 + Wi)ai 3Q(1 + Wi)&i 



(25) 



K 



n 



{VH - V)A t + Z>(1 + Wi) (k 2 a - 3r)) + P/i- 



C. Modified power spectra 



In the initial release of ISiTGR for both weak lensing tomography and ISW-galaxy cross-correlation measurements, 
the output 3D power spectrum Ps,s(k, z) was calculated by CAMB and then used to calculate the power spectra needed 
for lensing, P$ j $(fc, z) 7 and the power spectrum need for ISW-galaxy cross-correlations, P S2 ^{k,z), where $ = 
and A and B in Pa,b refer to the transfer functions used in the calculation of a given power spectrum. In this updated 
version of ISiTGR these power spectra are calculated directly in CAMB. A detailed description of the weak lensing and 
ISW-galaxy cross-correlation likelihoods can be found in [40j], but we will review the most relevant equations here. 

The lensing cross power spectrum between redshift bins k and I is given by 165] 



Jo a iX) \jk{X) 



'*•*(" 



,x , 



(26) 



with comoving radial distance x, comoving distance to the horizon x^, comoving angular diameter distance /if(x)j 
and where we have absorbed the usual extra terms into the P$,$. The weighted geometric lens- efficiency factor for 
the k th bin gk(x) given by 



9k(X) = / dxPk(X) — 7 / a ■ 



(27) 



corresponding to the galaxy redshift distributions p k - 
The ISW-galaxy cross power spectrum is written as [46 



a 



(£+l/2) 2 



(28) 



where b{z) is the galaxy bias, H(z) is the normalized selection function, and again we have absorbed the usual extra 
terms into the P, 
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IV. RESULTS 



A. Modified growth and cosmological parameters used 

For all results, in addition to the curvature parameter ttk and the MG parameters we vary the six core cosmological 
parameters: Qbh 2 and fl c h 2 , the baryon and cold dark matter physical density parameters, respectively; 9, the ratio 
of the sound horizon to the angular diameter distance of the surface of last scattering; r re i, the reionization optical 
depth; n s , the spectral index; and lnl0 lo A s , the amplitude of the primordial power spectrum. 

We use three different parametrizations of the MG parameters: two scale-dependent methods including a traditional 
binning method and a hybrid evolution method, both of which were discussed in our previous work [40(, and a scale 
independent method using a functional first introduced by [33]: 

1. For traditional binning both redshift, z, and scale, k, are binned in two bins creating a total of four bins. The 
z-bins are < z < 1 and 1 < z < 2 while for z > 2 GR is assumed valid. The fc-bins are simply k < 0.01 and 
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k > 0.01. This binning can be described functionally as 

Q(M) = i(l + 0, 1 (fc)) + ^(Q 22 (fc)-g 2l (fc))tanh 0nr> 



1 1 , , nN z-2 

-+2( 1 -^( fc )) tanh W 



(29) 



V(k,a) 



\{l + V Zl (k)) + \ {V z , (k) - V Z1 (k)) tanh 1^ 



with 



QM = 




HQi)- 


h ^° 2 


— Qi) tanh 


QzM = 


\{Q,- 


l- Qi) - 




— Q3) tanh 


v zi (k) = 








— T>ij tanh 


V Z2 (k) = 






"5^ 


— V3) tanh 



1 z — 2 

-(l - V z Jk)) tanh , 

2 V lK " 0.05 



fc-0.01 



0.001 ' 
fc-0.01 

0.001 ' 
fc - 0.01 

0.001 ' 
k - 0.01 

0.001 ' 



(30) 



(31) 



(32) 



This gives 8 MG parameters, T>i and Qi 



1, 2, 3, 4 which are varied. 



2. The redshift evolution hybrid evolution method is identical to that of traditionally binning ie Eq.(f2T)|). The 
evolution in scale, however, is described by a monotonic function written as [40], 



Qz x (k) = Qie-— +Q 2 (1 - e -— ), 

Qz 2 (k) = Q 3 e'^ +Qa(1 -e~om), 

V Zl (k) = V x e~^ +V 2 (1 -e"rar), 

V Z2 (k) = V 3 e~^ +D 4 (1 -e"^). 



(33) 
(34) 



This provides an exponential transition in scale between two parameter values in each redshift bin. Again a 
total of 8 MG parameters are varied, £>,; and Qi, i = 1, 2, 3, 4 



3. For the functional form method we vary the MG parameters Q and R according to [3c 

0(o) = (Q0-1K + 1, 

i?(a) = (i? -l)a s + l- 



(35) 
(36) 



The parameters Qo and Rq are the present day values of the MG parameters which go to 1 at early times. 
The parameter s describes the time dependence of the MG parameters and is marginalized over. To alleviate 
a known degeneracy in the parameter space between Qo and Ro, we vary Qo and Vq = Qa(l + Rq)/2 inferring 
Rq from the other two parameters. 



Correlation coefficients between Q,k and the MG parameters 




MG parameters evolved 


using traditional binning 




Qi 


Q 2 


Qs 


Qi 


Vi 








-0.1658 


0.0108 


-0.1114 


-0.1872 


-0.0752 


-0.3606 


-0.0327 


-0.0988 




MG parameters evolved using hybrid binning 


Qi 


Q 2 


Qs 


Qi 


V 1 


V 2 




Vi 


-0.0662 


-0.0840 


-0.1258 


0.1365 


0.0121 


-0.2480 


0.1025 


-0.2316 





MG parameters 


evolved using the functional form 


Qo 


V 


Ro 


-0.0166 


0.1428 


0.0400 



TABLE I: We list the correlation coefficients between ilk and the various MG parameters for the current observed data. We 
see that the MG parameters for all evolution methods are somewhat correlated with fl^- The functional form method shows 
overall the least amount of correlation between the MG parameters and the two cosmological parameters. The hybrid binning 
is the next least correlated of the methods, while parameters from the traditional binning method of evolution show the most 
amount of correlation with fife. 
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FIG. 1: We plot the 2D confidence contours for ilk and the MG parameters Qi and T>i, i = 1,2,3,4 from using traditional 
bins for k and z. As seen in Table U this evolution method has the most correlation between the MG parameters and ilk of the 
three evolution methods. 
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FIG. 2: We plot the 2D confidence contours for ilk and the MG parameters Qi and T>i, i = 1, 2, 3, 4 from the hybrid method to 
evolve the MG parameters. As seen in Table [I] this evolution method has a moderate amount of correlation between the MG 
parameters and ilk- 



B. Results from current data and correlations between ilk and the MG parameters 

Here we study the correlation coefficients between the MG parameters and il k . We choose to focus on this parameter 
because a quick examination of Eqs. §5§ and (TTU|) shows that the parameters should be correlated. Indeed, using Eq. 
(HOP in particular shows how this correlation is expected. The first term on the right-hand side of this equation has 
K, the spatial curvature, in the denominator while the MG parameter T> is in the numerator. An increased K and 
thus decreased Qk will (for a given range of wave number k) be balanced by a decreased T>. This is expected since an 
increased K (decreased f^) makes the denominator smaller and thus the whole first term larger. A decreased V will 
balance this effect. The same effect can be seen for the MG parameter Q if we rearrange Eq. © in a form similar to 
Eq. (fTU|) . The presence of this degenerate effect between the MG parameters and Slu is the reason we choose to focus 
on the said correlation in this paper. No other correlations are as obvious through an analytical examination. 

We first discuss correlation coefficients that were obtained from constraints using the latest cosmological data 
including; the WMAP7 temperature and polarization spectra (CMB) (46|, the matter power spectrum (MPK) from 
the Sloan Digital Sky Survey (SDSS) DR7 0, the ISW-galaxy cross-correlations [H,[49|], the refined HST-COSMOS 
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FIG. 3: We plot the 2D confidence contours for Qk and the MG parameters Qo and T>o, from using the functional form to 
evolve the MG parameters. As seen in Table [J this evolution method overall has the least amount of correlations between the 
MG parameters and the Qk of the three evolution methods. 



weak-lensing tomography I50l. WiggleZ BAO measurements [5l| , and the supernovae compilation of the Supernovae 
Cosmology Project (SCP) [52[ and references of other compiled supernovae (SN) therein. We also use the prior on 
H = 74.2 ± 3.6 km/s/Mpc given by [H, and a prior on the age of the Universe (AGE) 10 Gyrs<AGE< 20 Gyrs. The 
constraints on the MG parameters we obtain are all consistent with general relativity at the 95% confidence level. 

In Table U we list the correlations between fi^ and the MG parameters. We have used the standard definition for 
the correlation coefficient of two parameters, p x and p y : 

Corrip^) = (37) 

where the covariance of the two parameters is Cov(p x ,p y ) and their respective standard deviations are cr(p x ) and <t(p v ). 
We further illustrate these correlations by plotting the 2D confidence contours for the functional form, traditional 
binning, and hybrid evolution methods in Figs. UKU and|31 respectively. 

The existence of non-negligible correlations between £l k and the MG parameters indicates that curvature may 
possibly produce an apparent deviation from GR in the MG parameter space if the universe is assumed to be flat. 
To explore this possibility further, in the next section we generate future cosmological data with different curvature 
values and see how much the MG parameter constraints are affected by assuming a spatially flat background when, 
in fact, the background is curved. 



C. Constraints from simulated future data and the effect of assuming a spatially flat model on curved ones 

To explore the effect of assuming a flat background when the actual cosmology is curved (even slightly) we generate 
fiducial data with uncertainties achievable by future higher precision observations. It is worth noting here that the 
constraints obtained in this section are not meant to represent the constraining power of future observations, but 
rather to illustrate how the assumption of spatial flatness may affect constraints on MG parameters. The fiducial 
data sets are produced using a curved ACDM model so all of the MG parameters are set to their GR values during data 
generation. For each fiducial data set, we generated data for the following observations: CMB (Plank-like), matter 
power spectrum, weak-lensing, ISW-galaxy cross-correlations, Type la supernovae, and baryon acoustic oscillations. 
The CMB data was produced using CosmoMC and we used the included PYTHON script to add noise and make it 
compatible with the included Plank likelihood which was used during the fits. For ISW-galaxy cross-correlations, 
weak lensing and matter power spectrum data, the current galaxy distribution and window files were used. Random 
noise up to 5% of the true data value was then added to the data. For likelihood calculations uncertainties were 
reduced by a factor of 3 compared to current observations. BAO observations were given random noise of up to 
5% while supernova data points were allowed random noise up to 2.5 times current uncertainties. Uncertainties in 
these two data sets were not improved compared to currently available data, primarily because they probe only the 
expansion history of the universe and so are not important in constraining the MG parameters. For all fits, in addition 
to varying the MG parameters, the six core cosmological parameters are varied as well. A null test, where fifc is also 
varied, was performed on all fiducial data sets generated and the results of these fits in all cases recover the parameter 
values near those used to generate the fiducial data. This ensures that if indeed non-GR constraints are obtained for 
the MG parameters when we assume a flat universe, curvature is the cause and not a flaw in the production of the 
fiducial data. 

Two curved fiducial models were initially tested: one negatively curved model with = 0.01 and one positively 
curved model with VLf. — —0.02. These values were picked because they are the rounded WMAP7 95% confidence 
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FIG. 4: Plotted are the 68% and 95% 2D confidence contours for the MG parameters when a spatially flat model is used but 
the actual underlying universe has ilk = O.Of . TOP: Confidence contours for the MG parameters Qi and T>i, i = 1, ..4 when the 
traditional binning method is used. All parameter values are pulled to smaller values. MIDDLE: Confidence contours for the 
MG parameters Qi and T>i, i = 1, ..4 when the hybrid evolution method is used. Most of the parameter contours are pulled to 
smaller values. BOTTOM: Confidence contours for the MG parameters Qo and Do and Ro when the functional form binning 
method is used. The Qo — fo contour is pulled noticeably toward smaller parameter values. 



limits on fi^. The constraints on the MG parameters for these data sets is plotted in Figs. 0] and [3] respectively. 
Already for these values of Qk one can see that the GR point (shown where the two dot-dashed lines cross in each 
figure) is moving away from the best-fit point. In fact, in the case where Qk = —0.02, for the functional form evolution 
there is very nearly a deviation from GR at the 95% level. This shows that indeed ignoring curvature can possibly 
produce a false positive for a deviation from general relativity when higher precision data is used. 

To further explore how large an impact assuming the universe is fiat when in fact it is curved will have on constraints 
on the MG parameters four more fiducial data sets were included with ilk = ±0.05 and = ±0.1. Though these 
values for the curvature parameter ilk are well outside current constraints from all combined data, these constraints 
on ilk were obtained by using a ACDM model, and while in the ACDM cosmological model observations require 
a universe that is flat or very close to it, modified theories of gravity or inhomogeneous cosmological models may 
require the universe to be significantly more curved to fit observations. Additionally, given the fact that our simulated 
data sets do not represent the constraining power of future experiments, it is useful to use these larger values of 
the curvature parameter to illustrate how constraints from data sets with smaller uncertainties or more data will be 
affected. While the data sets we simulated and used may require these large curvature values for a significant apparent 
deviation to arise, future data sets with more precise measurements may not. For most of these values of ilk, we find 
for every parameter evolution method (traditional binning, hybrid method, or functional form) there is at least one 
MG parameter and in most cases more than one MG parameter that strongly deviates from its GR value. Almost all 
Q — V 2D parameter contours for these values of ilk show the GR, (1, 1), point outside their 95% confidence limits 
constraints. In Figs. [S]and[7]the 2D confidence contours for the ilk = 0.05 and ilk = —0.1 fiducial data sets are shown 
respectively. Interestingly the negatively curved model (ilk = 0.05) deviates from GR much more substantially than 
does the closed model for all evolution methods of the MG parameters. This could be due to the way the comoving 
angular diameter distance enters into both the weighting factor for the weak lensing as well as the way the wave 
number, k, is determined when calculating the lensing cross power spectrum. 
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FIG. 5: Plotted are the 68% and 95% 2D confidence contours for the MG parameters when a spatially flat model is used but 
the actual underlying universe has Qk = —0.02 TOP: Confidence contours for the MG parameters Qi and T>i, i = 1, ..4 when 
the traditional binning method is used. All parameter values are pulled to larger values. MIDDLE: Confidence contours for 
the MG parameters Qi and T>i, i = 1, ..4 when the hybrid evolution method is used. Most of the parameter contours are pulled 
to larger values. BOTTOM: Confidence contours for the MG parameters Qo and X>o and Ro when the functional form binning 
method is used. The Qo — 'Do contour is pulled noticeably toward larger parameter values. 
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MG parameters evolved using the functional form 



0.0289 



V 



0.0969 
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TABLE II: We list the correlation coefficients between Qk and the various MG parameters for a universe with fit = 0.05. 
In contrast to those obtained from current observational data, these correlation coefficients are almost all consistent with the 
trends observed in Fig. [8] when a flat universe is assumed but curved simulated future data is used. 



One other feature to notice in all four sets of plots is the directions which the constraint contours move in the 
parameter space. Looking at Eqs. © and (fTU)) . one would expect that for a universe with a positive value of f2 fc and 
thus a negative K value the assumption of a flat model would demand smaller values for the parameters Q and T> 
and vice versa. This is indeed the trend we see for the most part in Figs. IHE1EJ and [7] To further show how fiducial 
models with higher values of fi^ have lower best-fit MG parameters when a flat universe is assumed, in Fig. [5] we plot 
the best-fit MG parameter values versus flk for the various evolution methods. As was discussed above most of the 
parameters exhibit a negatively sloping trend. This seems to suggest that the observed behavior where most of the 
Q and T> parameters move the same direction in the parameter space could act as a signature of a false positive in 
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FIG. 6: Plotted are the 68% and 95% 2D confidence contours for the MG parameters when a spatially flat model is used but 
the actual underlying universe has Qk = 0.05. TOP: Confidence contours for the MG parameters Qi and T>i, i = 1, ..4 when 
the traditional binning method is used. All parameter values are pulled to smaller values and indicate a deviation from general 
relativity. MIDDLE: Confidence contours for the MG parameters Qi and T>i, i = 1, ..4 when the hybrid evolution method is 
used. Most of the parameter contours are pulled to smaller values and indicate a deviation from general relativity. BOTTOM: 
Confidence contours for the MG parameters Qo and Do and Ro when the functional form binning method is used. A strong 
deviation (due to the assumption of spatial flatness) from reneral Relativity is present in both contours. 



the future if flatness is assumed and a deviation from GR is detected. 

One can also notice that the trends discussed above do not match the behavior expected from the correlation coef- 
ficients obtained when using the current data. For the Q and T> parameters the behavior described above corresponds 
to positive correlation because by assuming flatness (for example, in the case of a positive fifc) we are forcing a lower 
value for VL^ and thus would expect a lower value for our MG parameters. We do not observe this behavior with the 
current data. This is most likely due to the large uncertainties in this data which allows a large number of models to 
fit quite well, and thus makes the calculation of the correlation coefficients less accurate. A more accurate description 
of the correlation coefficients can be obtained by using some that were calculated when running the "null" tests on 
our fiducial models. We show the correlation coefficients for the case when ilk = 0.05 in Table [TTl These correlation 
coefficients are almost all consistent with the trends we see in the constraints as well as the best-fits. 



V. CONCLUSION 



In this work we extended previous studies and the framework of modified growth (MG) parameters to test general 
relativity (GR) at cosmological scales in order to include spatial curvature in the models, as while current data when 
analyzed using the ACDM model points to a universe that is flat or very close to it, this constraint may not hold in 
modified theories of gravity or inhomogenous cosmological models. Using the latest cosmological data sets we explored 
the correlations between MG parameters and the curvature parameter J7fc finding that indeed there are non-negligible 
correlations. We next used future simulated data to explore whether assuming a spatially flat model on a spatially 
curved universe would affect the MG parameter constraints. We found that indeed, for our simulated data sets, such 
an assumption of flatness can cause tension with GR when using an f2fc as little as 0.02 away from the flat case 
and significant apparent deviations for \flk\ > 0.05 . Models with larger departures from flatness cause more MG 
parameters to deviate from GR with even larger discrepancies. We also found that negatively curved models deviated 
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FIG. 7: Plotted are the 68% and 95% 2D confidence contours for the MG parameters when a spatially flat model is used but 
the actual underlying universe has — —0.1 TOP: Confidence contours for the MG parameters Qi and T>i, i = 1, ..4 when 
the traditional binning method is used. All parameter values are pulled to larger values and indicate a deviation from general 
relativity. MIDDLE: Confidence contours for the MG parameters Qi and T>i, i = 1, ..4 when the hybrid evolution method is 
used. Most of the parameter contours are pulled to larger values and indicate a deviation from general relativity. BOTTOM: 
Confidence contours for the MG parameters Qo and X>o and Ro when the functional form binning method is used. The Qo — T>o 
shows a deviation from general relativity, while the Qo — Ro contour shows some tension with the GR point. 



more quickly and more significantly from GR when the assumption of flatness is made. 

As expected from the derived modified growth equations, for all approaches and for most of the bins, we find 
positive correlation coefficients between the MG parameters and 51^ from using simulated future data which is much 
more precise than what is currently available. The trends for the best-fit MG parameters versus fiducial values of flk 
used are found to be consistent with these correlations. Though the values of the correlation coefficients found are of 
course dependent upon the parametrization used, the signs of the correlations coefficients for the most part remain 
consistent across parametrizations. 

The results obtained in this analysis show that when using high-precision data from future experiments in order 
to test general relativity at cosmological scales or look for deviations from it, one must take into account the effect 
of curvature. This point may also be relevant when trying to test other, alternative theories of gravity. Indeed, 
our results indicate that the assumption of a spatially flat universe when performing these tests can bias the MG 
parameter constraints, leading to apparent deviations from general relativity. 
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